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Abstract. With the William Herschel Telescope in La Palma we made IJK S observations of an area of about 
40' x 30' of the Local Group galaxy Draco. This allows us to describe Draco's late-type stellar population across 
the whole galaxy at a photometric level 2 mag deeper than the 2MASS survey. We detected the red giant branch 
(RGB) and measured the magnitude of the tip of the RGB in the three bands. From that in the / band we obtain 
a distance modulus of (m — M)o = 19.49 ± 0.06(siat) ± 0.15(si/s), in excellent agreement with a measurement 
from RR Lyrae stars. The peak of the (J — K a )o histogram at different Mk b suggests that Draco has a mean 
[Fe/H]= — 1.95±1.26 while fiducial RGB tracks of Galactic globular clusters indicate a mean [Fe/H]= — 1.33±0.72 
where the error corresponds to the spread around the mean value. There are significant differences between the 
colour-magnitude diagrams of stars in the inner, medium and outer areas of the galaxy. A metal poor (Z= 0.0004) 
intermediate-age population (about 1.6 Gyr old) is clearly present and emerges in particular between 6' and 
12' from the centre of the galaxy. A few additional carbon star candidates have been identified from both their 
location in the colour-magnitude diagram and from an indication of variability. The large scale distribution of 
late-type stars is smooth but irregular in shape; this points at a variation of inclination with radius. 



1. Introduction 

The faint low surface brightness elliptical Local Group 
galaxy Dra co (a 2 oo = 17 h 20 m 18.1 s , <5 200 o = 57°55'13"; 
Piatek et al. 2002) is located at about 80±7 kpc (Aparicio 
et al. 120011 Nemec I1985fl from the Sun. Bonanos et 
al. (2004) derived a distance modulus of (m — M) = 
19.40 ± 0.02(stat)±0.15(syst). Its size is approximately 
35.5' x 24.5' and about 14 kpc along the line of sight 
(Aparicio et al. I2001f) . The foreground reddening (E(B — 
V) = 0.03 ± 0.01; Stetson fl979b|) is lower than expected 
for a galaxy at a galactic latitude b = 35°. The spread in 
magnitude among RR Lyrae stars within the galaxy shows 
no evidence for internal dust (Bonanos et al.[2004). 

The optical colour-magnitude diagrams show a well- 
defined horizontal branch (HB) and a red giant branch 
(RGB) without a signature of the RGB bump although 
deep enough observations have been performed (Bellazzini 
et al. I2002f) . The diagrams are similar to those of glob- 
ular clusters and indicate the dominance of an old and 
metal-poor population (Stetson lTW9a„ Baade & Swope 
119610 . The location and width of the RGB is consistent 
with a mean metallicity [Fe/H]= —2.0 ± 0.8 dex (Carney 
& Seitzer H986fl . Differences in temperature and spectral 
properties on the upper giant branch stars indicate the 
presence of both RGB and AGB stars (Stetson lTMTjl . The 
main-sequence turnoff corresponds to an age of about 18 
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Gyr and there is an almost negligible younger population 
(Stetson et al. H985J) . Observations with the Hubble Space 
Telescope (Grillmair et al. H998H showed that: (i) the stars 
formed primarily in a single epoch; (ii) the initial mass 
function is very similar to that in the solar neighbourhood; 
and (iii) there is a non-negligible fraction of blue stragglers 
(Carney & Seitzer, 1986). From simulations of the colour- 
magnitude diagram while accounting for the abundance 
pattern within the galaxy Ikuta & Arimoto l|2002|) con- 
cluded that over a long star forming period (> 3.9 — 6.5 
Gyr) the star formation rate was low (1 — 5% of that of 
the solar neighbourhood); the stripping of interstellar gas 
by the Milky Way Galaxy perhaps ended the star forma- 
tion in Draco. The intensity of diffuse H Q (Gallagher et al. 

2003) implies a mass of the ionised gas below 10% of the 
stellar mass. This limit is 10 times larger than the upper 
limit for the mass of neutral hydrogen (Young [1999). 

A few hundred variable stars have been identified 
in various surveys. The most numerous are RR Lyrae 
stars but there are also Cepheids (Baade & Swope IT9611 
Deupreee & Hodson [19771 Nemec 119851 Bonanos et al. 

2004) . The Sloan Digital Sky Survey covered the entire 
galaxy; cross-identifications with existing optical cata- 
logues gave 142 additional candidate variable stars (Rave 

et ai.mni). 

Star counts show that the Draco galaxy is elliptical 
with an eccentricity e = 0.29±0.02, a position angle of the 
major axis of 88° ± 3° and a cutoff radius in the long axis 
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around 26' ±2' fHodee lT^rjl Irwin & Hatzidimitriou lTM)! 
Odenkirchen et al. 1200111 . The Odenkirchen et al. paper 
shows that the spatial distributions of giants, horizontal 
branch, and sub-giant stars down to i= 21.7 within an area 
of 27° around the centre argues against tidal effects. This 
conclusion was also reached by Piatek et al. Q2002[l . Deeper 
optical photometry by Wilkinson et al. (2004) confirmed 
a break in the light profile at « 25' and a decline in the 
velocity dispersion outside 30' which imply the existence 
of a kinematically cold population in the outer part of the 
galaxy. On the other hand, there is no difference between 
the scale length of the distribution of old (> 9 Gyr) and 
intermediate-age (2 — 3 Gyr) stars and both populations 
were either formed under the same kinematic conditions, 
or any initial difference was afterwards erased (Aparicio 
et al-H!]}. 

Spectra of giant stars in Draco have been obtained 
since the early 1970s in order to derive the abundance 
of the heavier elements and other stellar parameters 
(Hartwick & McClure JTG7Q . The results showed that 
within the galaxy there is a spread in [Fe/H] of at least 
1.53 dex (Shetrone et al. 1200 lbl Shetrone et al. [TUM 
Zinn I1978f) and that giants in Draco are metal poor 
([Fe/Hjss —2.0) compared to giants in Milky Way globu- 
lar clusters ('Ganterna HUTBI Stetson HUES Wmnick USP - 
In fact the location of stars redward of the fiducial RGB 
sequence, indicate that these stars have a higher metallic- 
ity but still below [Fe/H]= —1.45 dex; they must be car- 
bon stars (Shetrone et al. 1200 lb"}) . Inhomogeneities in the 
[Ca/H] and [Mg/H] abundance also support an abundance 
spread (Winnick .2003] Lehnert et al. H992|l . Supernova ac- 
tivity may account for this spread and for the removal of 
about 90% of the original gas content (Smith I1984|l . A 
handful of carbon stars were identified by Aaronson et al. 
(1982) and by Azzopardi et al. 1|1986J) . One of those is a 
symbiotic binary (Munari 1991) and another is lithium- 
rich and probably a low-mass AGB star of low metallicity 
in its thermally pulsing phase after a few third dredge-up 
episodes (Dominguez et al. I2004f) . Olszewski et al. (|1995fl 
confirmed three carbon stars as velocity variables which 
implies that they are likely binary stars similar to CH stars 
in the Milky Way halo. The sample of confirmed carbon 
stars was increased to six by Shetrone et al. (20 01b|) . 

The mean velocity of the Draco galaxy, determined 
from different samples of stars, is — 293.3 ± 1.0 km/s with 
a dispersion of about 8.5 ± 0.7 km/s (Hargreaves et al. 
1996, Armandroff et al. [1995). This dispersion is signifi- 
cantly larger than the velocity dispersion of other galaxies 
in the Local Group, except for the Ursa Minor and Carina 
galaxies. The measured dispersion is not dominated by bi- 
nary stars (Olszewski et al. 1996) and the mass to light 
ratio (145i 71 ) is consistent with the presence of a large 
quantity of dark matter (Hargreaves et al.[l996). 

Our present work concerns the intermediate and old 
stellar population in Draco and uses new mosaic observa- 
tions in the /, J and K s bands. Observations and data 
reduction are described in Section 2 while in section 3 we 
discuss the distribution of sources in the colour-colour and 
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Fig. 1. Distribution of the mean-error as a function of 
magnitude (top) and colour (J — K s - middle and I — J 
- bottom) for sources in Table 1. Error bars indicate the 
standard deviation of the mean, they have been both ob- 
tained with a bin of 0.25 mag in IJK S and of 0.1 mag in 
colour. Lines connect the mean points of each bin: contin- 
uous (/ and colour), dotted (J) and dashed {K s ). 

colour-magnitude diagrams at different distances from the 
centre; we also discuss the contamination by foreground 
stars. Section 4 presents a determination of the distance 
to the galaxy, a discussion about variations in age and 
metallicity of the stellar population across the galaxy and 
a discussion of known carbon stars. Conclusions are given 
in Section 5. 

2. Observations & Data Reduction 

Observations in the /, J and K s bands were obtained 
with the 4.2 m William Herschel Telescope (WHT) on La 
Palma (Spain) during the same runs described in Cioni & 
Habing ( 2005 ) . Due to the large extent in the sky a mosaic 
of images was required in all three bands. In particular we 
obtained 11x9 images of 4.06' x 4.06' in the J and K s 
band and 6 images of 16.2'xl6.2' in the /-band. In the 
/-band we used exposures of 6s, 60s and 300s and J and 
K s images were obtained by co-adding 4 images out of 5 
of a single dithered position each with an exposure time of 
1.5s. In total we covered an area of about 40' x 30' centred 
at a = 17:20:12.4 and S = +57:54:55. 

The processing of raw images was done using the IRAF 
software and following the same procedure described in 
Cioni & Habing (2005J with a few minor differences: sky 
frames in the near-infrared (near-IR) wave bands were 
obtained averaging uncrowded target images; - the same 
fringe image was used to correct long and short /-band 
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exposures. Near-IR images were combined in 11 columns 
per band (9 images per column) while /-band images were 
combined in one single image using the spatial overlap. 
Each near-IR column was matched to the /-band image 
to find source counterparts. Note that there was no overlap 
among near-IR columns (= no overlap in right ascension), 
while the overlap in declination was quite large (1') to 
compensate for the effects of image distortions. 

Sources were extracted using the SExtractor program 
(Bertin & Arnout 1996). The flux corresponding to each 
source was calculated within an aperture of 5 pixels, a 
detection threshold of 2.0, an analysis threshold of 1.5. 
In the SExtractor configuration files we used these pa- 
rameters as well as the parameters that characterise the 
detector and default values for the remaining keywords. 
Sources were extracted first in each wave band separately. 
Afterwards J and K s detections were matched using an 
association radius of 6 pixels (6 pix x0.238"/pix w 1.4") 
to account for possible field distortions, although J and K s 
images were aligned before extracting individual sources. 
The nearest source was kept as a counterpart and the re- 
sulting catalogue was filtered to retain only sources with 
K$ < — (J — K s ) + 18.5 (this artificial cut is shown at 
faint magnitudes and red colours in Fig. [5] The cross- 
identification of /-band sources was made separately for 
each near-IR column. In practice an /-band column-image 
was extracted from the full /-band image to match a given 
near-IR column (i.e. by aligning common sources and 
trimming the area). Sources were matched using the same 
criteria but an association radius of 80 pixels. Different 
association radius were initially used and the resulting 
matches where checked against known sources. Finally 
wrong associations were removed by accepting only pairs 
for which 14 < / < 20 and -1 < (/- J) < 3. This criteria 
combines both position, magnitude and colour which guar- 
antees a statistical match of similar sources and accounts 
for field distortions. The /-band field of view is rather large 
(16.2' x 16.2') and distortions at the far edges are non- 
negligible. With a smaller matching radius we would have 
lost a few bright counterparts (those located at the edges 
of the field-of-view) . However because of the image align- 
ment and because we retained the closest match within a 
satisfactory colour and magnitude range good results were 
obtained. A few mis-matches, between the near-IR and I- 
band magnitudes might be present among the faint stars 
in the catalogue. For these sources the match between J 
and K s is more reliable compared to the match between 
JK S and /. 

The photometric calibration was performed for each 
observing night and it is described in detail by Cioni & 
Habing (2005). Photometric zero-points in the J and K s 
bands have been corrected for systematic shifts after cross- 
identification with 2MASS data. This process showed a 
systematic shift of 5.25" ± 1.75" between the astrometry 
of the two data sets were the error-bar indicates the disper- 
sion of the mean. Only sources within 1.75" of this value 
(248 sources) have been used to compare the photometry 
obtained within similar filters. We obtain: 



j = J2MASS + o.3o ± 0.05 
K s = K™ ASS + 0.30 ± 0.07 

These shifts (partly due to different filters) are smaller 
than those obtained by Cioni & Habing for the NGC 6822 
galaxy, probably because because Draco is not as crowded 
as NGC 6822, although both galaxies were observed dur- 
ing poor photometric conditions (i.e. thin cirrus). In the 
/-band we applied the same shift as derived for NGC 6822 
namely / = I BENIS - 0.1 ± 0.01. This is justified because 
the /-band observations of Draco and NGC 6822 were ob- 
tained during the same night under similar photometric 
conditions whereas our near-IR observations were spread 
over different nights. 

2.1. Catalogue 

Table H shows the first ten lines of the full table (ta- 
ble, dat) that is only available in the electronic edition 
of this article. The table contains 2570 stars detected in 
all three bands, /, J and K s : columns 1 and 2 list Right 
Ascension and Declination in degrees at the epoch J2000, 
columns 3, 4 and 5 list / magnitude, photometric error 
and SExtractor flag (1 = Bright neighbours or bad pixels 
affecting more than 10% of the integrated area, 2 = The 
object was originally blended with another), respectively; 
columns 6 — 8 and 9—11 contain the same information for 
the J and K s bands, respectively. The table lacks sources 
in the gap between pairs of /-band chips. We include only 
sources with SExtractor flag < 4 in all three wave bands. 
Figure shows the behaviour of photometric errors for all 
extracted sources as a function of magnitude in each wave 
band and as a function of / — J and J — K s colour. 

Crowding is unimportant in Draco and there is no issue 
about confusion. The catalogue is almost complete down 
to magnitudes 18, 16.8 and 17.2 in /,/ and K s but we 
may have missed a small percentage of sources because 
of missing adjacent tiles of the mosaic; this effect will be 
strongest in the /-band. These limits correspond to a de- 
crease in the number of sources towards faint magnitudes 
estimated from Fig. Sources detected in one single band 
or only in the two infrared bands are a very small percent- 
age of those presented in Table 1; the catalogue includes 
all stars in the upper RGB magnitudes (Sect. 3). These 
single band sources or those detected only in / and K will 
be made available upon request to the first author. 

3. Results 

3.1. Magnitudes and colours 
3.1.1. /, J and K s histograms 

Figure [21 shows the histogram of the number of sources 
detected in all three wave bands as a function of magni- 
tude. The discontinuity produced by the TRGB is clearly 
visible in the /-band but it is not obvious in the J and 
K s magnitude distributions. Other discontinuities may be 
present; they will be considered in some detail in Sect. 3.3. 
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Table 1. Draco catalogue of sources detected in three wave bands. 



a S I aj fj J a j fj K s a Kg i Ks 

260.089081 57.691799 18.04 0.13 17.05 0.17 16.16 0.07 

260.090729 57.697617 19.07 0.21 17.05 0.17 16.16 0.07 

260.048615 57.717018 16.65 0.07 16.06 0.11 15.37 0.06 

260.084015 57.717052 17.99 0.13 16.90 0.16 16.08 0.07 

260.066681 57.717663 18.63 0.18 18.40 0.27 18.09 0.08 

260.086700 57.718666 17.96 0.13 16.90 0.16 16.08 0.07 

260.045319 57.720440 17.92 0.13 16.06 0.11 15.37 0.06 

260.008881 57.723244 18.98 0.21 18.38 0.27 17.14 0.08 

260.003357 57.725845 17.80 0.12 18.38 0.27 17.14 0.08 

260.065857 57.726101 18.55 0.17 18.40 0.27 18.09 0.08 
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Fig. 2. Histogram of the number of sources matched 
among the three wave bands as a function of magnitude: 
I (continuous line), J (dotted line) and K s (dashed line). 
Each bin is 0.5 mag. 

The precise location of the TRGB, derived using the same 
technique as in Cioni et al. (|2000bi) . is at I = 15. 54 ±0.06, 
J = 14.89±0.10and K s = 14.07±0.08. The algorithm that 
determines the TRGB position starts with an estimated 
magnitude and a magnitude range. In the I band we used 
I = 15.5 as from the discontinuity in the magnitude distri- 
bution and a range of 0.5 mag. This is consistent with the 
value derived by Bellazzini et al. (2002) and accounts for 
possible variations between the two photometric systems. 
The TRGB in the J and K s bands has not been derived 
prior to this work, and so we had to derive the estimate 
ourselves. Using the isochrones by Girardi et al. (2000) 
for a 9 Gyr old metal poor population [Z = 0.0004; this 
corresponds to about [Fe/H]= —1-7) we derived the dif- 
ference between the TRGB magnitude in the /-band and 



those in the J and K s bands; they are 0.8 and 1.5 mag, 
respectively. This results in initial values for the TRGB of 
J = 14.7 and K s = 14.0 and the same range of 0.5 mag. 
Madore & Freedman (1995) postulated that at least 100 
stars in the upper 1-magnitude bin of the RGB are needed 
for a robust determination of the TRGB magnitude (see 
also Fig. 10 in Bellazzini et al. I2002f) and we do fulfil this 
criterium: we dispose of more than 200 stars in the up- 
per 1-magnitude bin in the / and J bands and more than 
100 in the K s band. Errors associated to each TRGB de- 
termination are only statistical errors. Systematic errors 
due to errors in the zero-point and in the calibration of 
the magnitude scale amount to 07 = 0.09, <j,j = 0.11 and 
<jk s — 0.15. We have not removed potentially foreground 
stars; for Draco this is a difficult task and foreground stars 
are expected to have a continuous distribution in magni- 
tude and thus have no effect on the location of the discon- 
tinuity that is the TRGB. Further considerations about 
foreground stars are given in Sect. 3.4. 

3.1.2. (/ — J, /) colour-magnitude diagram 

The colour magnitude diagram (/ — J, /) is shown in Fig. 
EJ It is a sparse distribution of sources and character- 
istic features of such diagrams cannot easily be traced. 
The location of the TRGB marks the distinction between 
fainter sources that belong either to Draco or to the fore- 
ground and brighter sources that most probably belong 
to the foreground except those with I — J rs 1 that may 
be AGB stars in Draco. Note that sources with I w 13 
approach the saturation limit of our observations. Only 
sources with 07 < 0.18 have been plotted, this corresponds 
to sources with I < 18.5 that represent also the most re- 
liable matches among the three photometric bands dis- 
cussed in this paper. 

3.1.3. (I — J, J — K s ) colour-colour diagram 

The distribution of sources in the colour-colour diagram 
(Fig. shows no clear-cut separation between sources of 
a different type as was the case for the Magellanic Clouds 
(Cioni et al. EZHMjl and for NGC 6822 (Cioni & Habing 
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Fig. 3. Colour-magnitude diagram of stars detected in 
three wave bands with ui < 0.18 mag. The location of 
the TRGB derived in this work is indicated. 

2005} . As in Fig. [3] we have plotted only sources with 
a i < 0.18. Most of the sources are concentrated at about 
I — J = 0.5 and J — K s = 0.75. There are usually very 
few foreground stars with J — K s > 1 and I — J > 1, 
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Fig. 4. Colour-colour diagram of stars detected in three 
wave bands with aj < 0.18 mag. Regions of expected 
RGB, C-type and M-type AGB stars as well as Galactic 
foreground stars are indicated. Known carbon stars are 
marked with thick dots. 



although the expected region of RGB stars clearly over- 
laps with that of foreground stars. Only two known carbon 
stars seem to fall within the expected C-type region, some 
are much fainter while one is considerably brighter and 
approaches perhaps the region where M-type AGB stars 
could be located. It is important to remember that the 
colours of stars of a different type, and therefore the ex- 
pected regions in the colour-colour diagram, depend on 
the age and on the metallicity of the stellar population. 

3.1.4. Near-IR colour-magnitude diagram: confirmed 
members and non-members of Draco 

The colour-magnitude diagram (J—K s , K s ) of Draco stars 
(Fig. [SJ shows a well defined and broad red giant branch. 
To identify the main stellar populations in this diagram 
we cross-correlated our catalogue (Table 1) with confirmed 
carbon and RGB star samples and with non-member 
sources as given by Shetrone et al. I|2001a|) . We also cross- 
correlated our catalogue with foreground dwarfs identified 
by Olszewski et al. 111)9.5) and with the catalogue of con- 
firmed members and non-members by Armandroff et al. 
( 1995i). These cross-identifications have been made using 
only coordinates, except for the known carbon stars that 
have also been checked in the corresponding images, and 
a few mismatches might have occurred, especially among 
the faintest counterparts. In Fig.|S]we used colours to mark 
stars identified by agreement in position and stars that 
are saturated in the /-band (possibly foreground stars); 
the colours are better visible in the on-line version of the 
article than in its printed version. The sensitivity limit of 
previous observations limited our identifications to sources 
with K s < 15. Among sources brighter than the TRGB 
the brightest magnitudes are likely foreground stars. Many 
of these have saturated the /-band photometry and have 
predominantly (J — K s ) < 0.8. Redder sources could be 
intermediate-age AGB stars. One likely thermally pulsat- 
ing carbon star has been identified in this region (Sect. 
4.2). Along the main RGB sequence confirmed members 
and non-members are equally distributed and do not sin- 
gle out a specific region in the diagram. Figure shows 
also the histogram of the number of sources as a function 
of colour for the full sample (Table 1), for the confirmed 
foreground stars and for the sample that remains after 
subtraction of the foreground. The foreground component 
dominates at colours (J — K s ) < 0.8 and becomes negligi- 
ble at redder colours. Once the foreground stars have been 
removed from the full sample the distribution of sources 
looks fairly symmetric and can be described by a Gaussian 
with /i = 0.74, a = 0.24. 

3.2. Foreground stars from 2 MASS 

Despite the availability of the precise location of some 
foreground stars discussed above (Armandroff et al. 1995 
Olszewski et al. I199fil Shetrone et al. I2001a[) . they can- 
not be used to estimate the overall foreground contamina- 
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Fig. 5. (Left) Near-IR colour-magnitude diagram of stars in Draco. Known carbon stars, known RGB stars (Shetrone 
et al. I2001a|) and stars that are confirmed not to be members of the galaxy (Olszewski et al. IT9951 Armandroff et al. 
1995) arc indicated with different colours. (Right) Histogram of the (J — K s ) colour for sources with K s < 15.5 (thin 
continuous line), confirmed non- member sources (dashed line) and the subtracted histogram (thick continuous line). 
The best fit Gaussian has fx = 0.74, a = 0.24 and FWHM= 0.56. The width of each bin is 0.1 mag. 



tion because their parent sample was not complete. Infact 
Olszewski et al. (1996) observed the brightest and reddest 
objects selected from the literature including additional 
giants near the empirical RGB of Draco. 

We extracted from the 2MASS catalogue sources 
within an area of 0.2 deg 2 0°.4 to the North of the cen- 
tre of Draco. This region is sufficiently close to Draco 
to provide a comparable sample of foreground stars but 
far enough to contain very few, if any, genuine Draco 
sources. Absolute magnitudes have been assigned suppos- 
ing that they have the same distance modulus as the stars 
in Draco and that the absorption is equal to that dis- 
cussed in Sect. 4.1. Figure H3 shows the distribution of 
the stars in the off-set areas in magnitude, colour and 
in the colour-magnitude diagram; the off-set areas were 
of two sizes that coincide with the IN region and the 
MED region discussed in the next Sect. 2MASS obser- 
vations are not sensitive enough to a foreground contribu- 
tion directly comparable with our observations but they 
indicate how the foreground stars are distributed on the 
sky. Foreground stars cover almost homogeneously a re- 
gion with 0.2 < ( J — K s )o < 0.9 independently of their 
Mk b mag. Their number increases in proportion to the 
surface area. Sources with (J — K s ) > 0.9 are close to 
the detection limit of 2MASS and have therefore larger 
errors. A comparison between Fig. [2] and Fig. [B] shows no 
gap at Mk b ~ —5 (K s « 14); this confirms our earlier 
assumption that one may find the TRGB without first 
removing the foreground. The gap visible in the colour 



histogram ((J — K s ) w 0.55) may have been caused by 
differences among the nature of foreground stars (Cioni et 
al. l2000al Nikolaev & Weinberg EES} , see also Fig.0 

3.3. Distribution in three concentric areas 

We take for the centre of the galaxy the location de- 
rived by Piatek et al. H2002J1 . We have converted Right 
Ascension and Declination into circular coordinates and 
have divided the area surveyed into three regions: p < 0.1 
(inner=IN), 0.1 < p < 0.2 (medium=MED) and p > 0.2 
(outer= OUT); see Fig. Sources are distributed fairly 
homogeneously across the galaxy and do not cluster in 
the centre. However, there is a non-negligible difference 
between the near-IR colour-magnitude diagram of stars 
in three aspects: (i) a variation in the expected location of 
the TRGB, (ii) a variation in the number of stars brighter 
than the TRGB and (iii) a variation in the colour of the 
RGB. 

The discontinuity attributed to the TRGB from 
the distribution of all sources corresponds to Mk s — 
—5.42 (using the distance modulus derived in Sect. 4.1). 
Although the number of objects around this magnitude 
is smaller within each region it is possible to distinguish 
other discontinuities. These are indicated with horizontal 
lines in Fig. [7] Considering for example the discontinuities 
at M Ks = -4.25 (IN), -4.75 (MED) and -5.25 (OUT) 
one may conclude that there is a variation of about 1.0 
mag between the inner and the outer region in the sense: 
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Fig. 7. (Upper) Distribution of Draco sources in the plane of the sky. Right ascension and declination have been 
converted into x and y coordinates with the centre of the galaxy (Piatek et al. (2002) as origin. We consider three 
different regions (inner (IN p < 0.1), medium (MED, 0.1 < p < 0.2) and outer (OUT, p > 0.2). The three histograms 
of absolute magnitude Mk s and the three of colour (J — K s ) are shown with filled, shaded and empty histogram bins, 
respectively. The width of each bin is 0.4 mag in Mk s and 0.15 mag in (J—K s ). (Lower) The near-IR colour-magnitude 
diagram of sources in the same three regions, IN, MED and OUT. Horizontal lines indicate the approximate position 
of discontinuities in the star counts. 



the TRGB in the inner region is fainter than the TRGB in 
the outer region by about 1.0 mag. The photometric error 
of each single source is dominated by the systematic error 
in the distance modulus, but as we are interested in varia- 
tions in magnitude the only important error is that in the 
apparent magnitude; this is about 0.07 mag (cfg. Fig.^J. 



If, in the worst-case scenario, we assume that all sources in 
the OUT region of Fig. [7| belong to the foreground, then 
we find that about 15% and 45% of the sources in the 
IN and MED region should also belong to the foreground; 
most of them would have Mk b > — 4. 
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Fig. 6. (Upper) Magnitude and colour distribution of fore- 
ground sources extracted from the 2MASS catalogue in 
the vicinity of Draco. Filled and shaded histograms rep- 
resent number counts in an an off-set area equivalent to 
the IN and MED area in Fig. respectively. The width 
of each bin is 0.4 mag in Mr- 8 and 0.15 mag in J — K s . 
(Lower) The near-IR colour-magnitude diagram of fore- 
ground sources covering the same areas. 

The histogram of the (J — K s )q shows that the peak 
increases by about 0.1 mag by going from IN to MED to 
OUT. This is about equal to the error in the colour of 
each single source, but the error on the peak of each his- 
togram decreases as 1/y/N and is of the order of 0.01— 0.02 
mag. Unless the foreground distribution strongly alters the 
shape of these histograms (cfg. Fig. |SJ) we may consider 
the increase of 0.1 mag as real. 

In the colour-magnitude diagram of the IN region the 
RGB clearly shows a jump in the number of sources 
around Mk b = —4.25, although the diagram clearly con- 
tains brighter magnitudes at approximately the same 
colour until Mk s = —6. There are a few, even brighter 
sources but they have bluer colours and were found to be 
the foreground stars discussed below. In the MED region 
a major discontinuity occurs around Mk„ = —4.75 as well 
as a cut off in the number of sources around Mk s = — 6. 
Sources with —4.75 > Mk b > —6 trace a clear branch 
about 0.4 mag wide and have bluer colours compared to 
the main RGB traced by fainter sources. Their colour is 
also bluer than stars with similar magnitudes in the IN re- 
gion. These sources could be metal poor intermediate-age 
AGB stars. Some or all of the very bright sources could be 
foreground sources. In the OUT region there is a major 
discontinuity around Mk s = —5.25 and perhaps a fainter 
discontinuity around Mk s — —4.5. At much brighter mag- 
nitudes (Mk b < —6) sources are distributed rather homo- 




260.000 
Right Ascension 



Fig. 8. Logarithmic and smoothed density distribution of 
stars detected in Draco in 32 x 22 bins of about 1'.44 
square each. Darker regions correspond to a higher density. 
Contours are from 3 to 11 in steps of 1. 

geneously within a colour range of about 1 mag. They 
most probably belong to the foreground, although a few 
AGB stars may hide among them. 

3.4. Spatial distribution 

Figure [S] shows the number density of sources in table.dat 
in 32 x 22 bins where a single element corresponds to about 
1'.44 square, which is equal to about 27 pc. The map cov- 
ers an area of about 40' x 30' and this is smaller than 
the cut-off radius of the density profile at optical wave- 
lengths (25' - 28') (Hodge EHl Irwin & Hatzidimitriou 
ITMBI Odenkirchen et al-HH Wilkinson et al. 127701)) . The 
source density in Fig.[S]has been smoothed using a box car 
function of width = 2 prior to the construction of the grey 
scale image where higher concentrations of sources are in- 
dicated by darker regions. The distribution, strongly dom- 
inated by RGB stars, appears regular and elliptical with 
an ellipticity e = 0°.27 (where e = 1 — b/a; a = 13'. 2 is the 
semi-major-axis and b = 9'. 6 is the semi-minor- axis); and 
PA w 90°. Intermediate contours (at 5, 6 and 7) show an 
enhancement of the number of sources towards the NW 
and the SE corners, suggesting a variation in the inclina- 
tion of the elliptical structure with radius. 

4. Discussion 

4.1. Determination of the distance to the Galaxy 

The /-band magnitude of the TRGB depends only weakly 
on age and metallicity (Salaris & Cassisi 1997) and can 
thus be used as a standard candle; see Lee et al. I|1993ll . We 
adopt the interstellar extinction as measured by Stetson 
H1979b|> of E(B - V) = 0.03, and the absorption in the 
/, J and K s band derived using Cardelli's (1989) law. 
These are: Aj = 0.05, Aj = 0.02 and A K = 0.01. 
Combining the apparent /-band TRGB magnitude de- 
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rived in the previous section with the most recent abso- 
lute calibration of the TRGB magnitude in the /-band 
by Bellazzini et al. (j2UU2"jl . M/ = -4.04 ± 0.12, we obtain 
(m-M)o = 19.49±0.06(stat)±0.15(sys); this corresponds 
to a distance of 79 ± 10 kpc. The systematic error quoted 
includes the error on the absolute magnitude and on the 
zero-point in the /-band (ctj = 0.09; Cioni k. Habing[2005 ) 
as well as the error from the calibration of the photometry 
scale {a i = 0.02). 

Our determination of the distance modulus is in very 
good agreement with the most recent determination of 
the distance to the galaxy by Bonanos et al. I|2004|l from 
observations of RR Lyrae stars. It is about 0.3 magni- 
tudes smaller than the value obtained by Bellazzini et 
al. 1)2002(1 who already noticed the discrepancy between 
their result and other measurements in the literature 
(e.g. Mateo IT998P and they attributed it to their adopted 
zero-age horizontal branch distance scale (Ferraro et al. 
1999). Bellazzini et al. determined for the TRGB the value 
/ = 15.8 ± 0.2 (their Fig. 12), and this is 0.26 mag larger 
than ours, an acceptable difference when all uncertainties 
are considered. It may be that a shift in the photomet- 
ric scale is responsible for this discrepancy but on the 
other hand Bellazzini et al. I|2002|l did not have enough 
sources to safely determine the TRGB position although 
they carefully introduced artificial stars. Our results of 
the distance to the galaxy agree with the determination 
by Aparicio et al. I|2001[l and by Nemec l|1985f) based on 
B and R band photometry and observations of RR Lyrae 
stars, respectively. 

4.2. Carbon stars and other AGB stars 

Shetrone et al. (f2QDTa|) spectroscopically observed a sam- 
ple of stars located on the red side of the fiducial RGB 
in the {B — V, V) colour-magnitude diagram. They in- 
creased the number and membership of known carbon 
stars to 6 and identified non-members by measuring ra- 
dial velocities. The photometry from our study and from 
the 2MASS counterparts are indicated in Table El In this 
case the cross-identification between our catalogue and the 
2MASS catalogue has been checked on the observed im- 
ages. Source names are those commonly used except for 
C5 and C6 that have been introduced by us for complete- 
ness. Four out of these six carbon stars are located around 
the TRGB. The detection of Li in the brightest of them 
(Dominiguez et al. I2004JI suggests it is a thermally pulsat- 
ing AGB star. Figure [3] shows that around the location of 
this star there is a handful of stars too red to be foreground 
stars (J — K a > 0.8), that could eventually also be AGB 
stars. Spectra are needed to confirm their nature. If the 
population of Draco has a large spread in metallicity the 
TRGB will either be smeared out or be found at different 
K s magnitudes in different locations; it is not necessarily 
true that the fainter carbon stars will be below the TRGB 
(Sect. 3.3). It is difficult to distinguish oxygen-rich AGB 
stars, if any is present, only from photometric observations 



of a galaxy with a low AGB component. Oxygen-rich AGB 
stars will be located above the TRGB and at J — K s ss 0.8 
where foreground stars are also present or at J—K s > 0.8, 
where carbon stars might be present. There are too few 
sources in this region in Fig. below we give an alterna- 
tive to find AGB candidates. 

AGB stars often experience variations in brightness 
with a long period and a large amplitude. With only two 
infrared epochs we cannot determine the periodicity of 
a magnitude variation but we can obtain an indication of 
the minimum amplitude that such a variation would have. 
Figure [5] shows the difference between J and K s magni- 
tudes from the observations described in this paper and 
from the measurements extracted from the 2MASS cata- 
logue. The photometry given in Table has been already 
adjusted to the 2MASS photometry scale. For a signifi- 
cant number of sources A J > 0.2 mag and AK S > 0.2), 
these are candidate variable stars and their position in 
the near-IR colour-magnitude diagram is indicated in the 
right panel of the figure. About half of the sources differing 
in J-magnitude between the two catalogues are brighter 
than the TRGB, up to 2 mag, while the other half extends 
to 1 mag below the TRGB. These sources with magnitudes 
around the TRGB also have red colours. They are likely 
AGB stars and suitable targets for a spectroscopic study. 

4.3. Variations in age and metallicity 

4.3.1. The population in the whole area 

Dereddening the near-IR magnitudes and applying the 
distance modulus obtained in Sect. 4.1 it is possible 
to show that the mean metallicity of Draco is [Fe/H] = 
— 1.95 ± 1.26 (where the error is the variance of the mean). 
Figure ED shows the distribution of Draco sources in the 
((J — K s )q, MkJ colour-magnitude diagram. These are 
the same sources shown in Fig. and foreground stars 
are still included, which as can be easily judged from 
their location, do not affect the following discussion (at 
least at magnitudes where they have been confirmed). 
Horizontal lines indicate the location of photometric in- 
dices that may be used to determine [Fe/H] (Valenti et 
al. 1^004). In Fig.^|we show the histogram of all sources 
fainter than the TRGB (RGB stars only) and three his- 
tograms of stars within 0.05 mag from the reference mag- 
nitudes Mk s — —5, —4 and —3. The first histogram shows 
a fairly symmetric distribution centred at (J — K s )q = 0.6 
with FWHM= 0.2. The other three histograms peak at 
(./ — K s ) = 0.4, 0.8 and 0.6, respectively, and are not 
symmetric. According to Valenti et al. I|2004|l the peak 
colours correspond to [Fe/H]= —3.3, —0.8 and —1.76, for 
the three histograms, respectively. These results do not 
allow us to safely determine the slope of the RGB, espe- 
cially considering that our observations, though about 2 
mag deeper than 2MASS observations, do not reach the 
base of the RGB (about 7-8 mag below the TRGB). 
They confirm, however, that a metallicity spread within 
the galaxy does exist that amounts to about 1.26 dex (vari- 
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Fig. 9. (Left): Magnitude difference, J (lower) and K s (upper), between the near-IR photometry presented in this 
article and the 2MASS photometry for the same sources in Draco. Variable sources depart from a line at zero values. 
Note that that WHT photometry has been adjusted to 2MASS photometry before calculating the magnitude difference 
as in Sect. 2. (Right): Near-IR colour-magnitude diagram of Draco (WHT sources) where thick dots indicate stars 
with more then 0.2 difference in both the J and K s bands. 



Table 2. Near-IR photometry of known carbon stars 



Source 


I 




J 


CTJ 


K s 


<JK S 


j2MASS 


C J2MASS 


jj2MASS 


<Jjj2MASS 


ry-2MASS 
^ s 


(7 K 1M ASS 


2MASS - ID 


461 


16.31 


0.06 


14.58 


0.06 


13.71 


0.04 


14.77 


0.05 


14.15 


0.05 


14.06 


0.06 


17194237 + 5758376 


68 


16.81 


0.08 


16.26 


0.12 


15.69 


0.08 


16.35 


0.10 


15.65 


0.12 


15.56 


0.19 


17195726 + 5755042 


3203 


15.78 


0.05 


14.76 


0.06 


13.90 


0.05 


14.38 


0.03 


13.71 


0.04 


13.46 


0.04 


17195764 + 5750054 


J 


17.75 


0.12 


17.41 


0.20 


17.35 


0.08 


14.92 


0.04 


14.22 


0.05 


14.05 


0.06 


17200068 + 5753464 


3237 


15.81 


0.05 


15.42 


0.08 


14.60 


0.06 


15.29 


0.06 


14.54 


0.06 


14.39 


0.10 


17203352 + 5750196 


578 


16.49 


0.07 


16.14 


0.12 


15.54 


0.07 


15.90 


0.10 


15.49 


0.14 


15.66 


0.29 


17203883 + 5759346 



ance of the mean) around a mean value of — 1.95 dex. Both 
the mean and the spread agree with other measurements 
in the literature (e.g. Hartwick & McClure 1996;, Carney 
& SeitzerHHEEl Shetrone et al. l2001bl ZinnQHEI Winnick 

However three facts need to be considered: (i) there 
is a systematic error in the (J — K s )q colour due to the 
adjustment of the photometry to the 2MASS catalogue, 
(ii) metal poor AGB stars are present at —5 > Mjf s > —6 
and (J — K s )o w 0.6 (Sect. 4.3.2) and (iii) the foreground 
contribution at faint magnitudes (Mk s < —4) has not 
been determined. Though the statistical error in the colour 
is about 0.2 mag, because of the larger error in the J band 
compared to the K s band, it decreases with the number 
of sources in each bin as 1/y/N and it becomes on average 
0.05. The systematic error on the colour amounts to about 
0.09 mag which corresponds to an average spread in the 



calculation of [Fe/H] using the Valenti et al. I|2004|l indexes 
of about 0.57 dex. If indeed sources with —5 > Mrt s > —6 
and (J— K a )o ~ 0.6 are metal poor intermediate-age AGB 
stars they should be excluded from the analysis of the 
slope of the RGB. This is equivalent to reduce the peak 
in the histogram of sources with Mk s = —5. We have 
seen in Sect. 3 that Galactic foreground with Mk s > —4 
stars have (0.2 < J — K a < 0.9) colours; their number is 
expected to increase at fainter magnitudes and their colour 
to cover a larger region if photometric errors become also 
large (which is usually the case approaching the sensitivity 
limit of the observations). Therefore we conclude that the 
colour of the RGB both at Mk b — —4 and Mr, = —3 are 
uncertain. 

Another way to estimate the metallicity from the 
( J — K s , MjfJ colour-magnitude diagram is to compare 
it with fiducial lines of Galactic globular clusters with 
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Fig. 10. (Left) Colour-magnitude diagram ((J — K s )o, Mr-J of stars in Draco. Thick slanted lines indicate the fiducial 
RGBs, from Ferraro et al. (2000 ), of three globular clusters (from left to right: M15, M4 and NGC 104). Their average 
metallicity is also indicated. Thick horizontal lines indicate reference magnitudes around which histograms in colour 
have been constructed with a bin size of 0.2 mag (right). These show the distribution of sources below the TRGB 
(Mk b — —5.43), and within ±0.05 mag of Mk b — —5, —4 and —3, respectively. The binning is comparable to the size 
of the error in the colour; it decreases with 1/y/N in the histogram. 



known metallicities distance and reddening (e.g. Ferraro 
et al. l2000[) . This comparison with clusters M15, M4 and 
NGC 104 is also shown in Fig. ^| These clusters have 
[Fe/H]= -2.12, -1.19 and -0.7, respectively. The thick- 
ness of the rather dispersed diagram of Draco sources is 
fitted well by the RGB of these clusters indicating an av- 
erage [Fe/H]= -1.34±0.72 dex. A fit of the TRGB is also 
clear indicating that the majority of Draco stars are as old 
as Galactic clusters. 

4.3.2. The population in three concentric areas 

Theory predicts that in two samples of equal age the 
TRGB is brighter in the sample with higher [Fe/H]. If 
the metallicity is the same, the TRGB is fainter in the 
older sample. Using the isochrones by Girardi et al. ( 2000) 
we have calculated the difference in the TRGB position 
between populations of 3 and 18 Gyr old and between 
populations with Z= 0.0004 ([Fe/H]= -1.7) and Z=0.001 
([Fe/H]= —1.3); these values were selected because we ex- 
pect Draco to be dominated by an old and metal-poor 
population. The difference in TRGB magnitude equals 
0.244 ± 0.003 mag when the two samples differ in age and 
not in metallicity and to 0.284 ±0.003 mag when they dif- 
fer in metallicity but not in age. On average a difference of 
about 0.26 ±0.3 is expected by a variation of either metal- 
licity or age. Between a pair of regions (e.g. IN-MED) we 
measure a difference of about 0.5 mag in the location of 



the TRGB (Fig. [TJ and this can only be explained by both 
a variation in age and metallicity. The colour of the RGB 
is also affected by variations in age and metallicity: at a 
constant metallicity younger stars have bluer RGBs while 
at a constant age the RGB of metal rich stars is redder 
than the RGB of metal poor stars. The average colour 
variation in either age or metallicity expected from the 
same populations discussed above amounts to 0.10 ± 0.03 
mag. This is of the same order as the variation in colour 
between a pair of regions and can thus be explained by a 
variation in age or metallicity only or by a combination 
of both. However, the difference in colour between the IN 
and the OUT regions is twice larger and must be due to 
both age and metallicity. Another non-negligible contri- 
bution to variations of the TRGB location is of course 
a variation in distance within the galaxy. Aparicio et al. 
(2001) suggested that Draco extends about 14 kpc along 
the line of sight. This corresponds to a variation of about 
0.8 mag in the TRGB position! 

We now discuss the apparent discontinuities in the star 
counts shown in the lower part of Fig. [7| The brighter 
discontinuity in the IN region is equally well described 
by Z= 0.0004 and age = 18 Gyr as by Z= 0.001 and 
age of 3 — 18 Gyr while the fainter discontinuity, if due 
to RGB stars at the TRGB, is better represented by a 
slightly more metal-rich population of comparable age. It 
is important, however, to be very careful in deciding which 
isochrone best fits the distribution of sources because pho- 
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tometric errors in Mk s are dominated by the systematic 
error in the distance modulus (0.15 mag) while errors in 
the (J — K s )q are dominated by the error in the J magni- 
tude and affected by a systematic error of about 0.09 mag 
(Sect. 2). These errors do not allow us to make firm state- 
ments about an absolute age or metallicity in a given re- 
gion but they permit the analysis of differences among the 
regions. In the MED region stars that produce the bright 
discontinuity are bluer than stars with similar magnitudes 
in the IN region. This means that in the MED region they 
are either AGB stars with Z= 0.0004 and age = 1.6 Gry 
or RGB stars with Z= 0.001 and age = 0.1 Gyr. There 
is no evidence of such a young population of giant stars 
from optical studies although a burst of star formation 
may have happened 2 — 3 Gyr ago (Aparicio et al.|2001). 
The fainter discontinuity can be explained by a population 
with Z— 0.001 and age = 18 Gyr while a population with 
Z= 0.0004 is bluer than the observed RGB unless a sys- 
tematic photometric shift is needed, i.e. the difference is 
due to measurement errors. The two discontinuities in the 
OUT region, because of the progressive reddening of the 
RGB, could be due to a population of increased metallic- 
ity (that could be as young as a few Gyr). The thickness of 
the RGB around Mk b — —4 suggests the presence of also a 
metal poor population, but at these faint magnitudes and 
blue colours foreground stars have not been studied and 
their distribution and number may bias this conclusion. 

Summarising: if most of the Draco population is as old 
as the Universe (13.7 ±0.2 Gyr; Bennett et al. 2003) then 
its metallicity increases with radius, which is against the 
usual trend. Otherwise if the spread in metallicity is not 
very large then the variation in age must be large and must 
increase with radius in agreement with the usual trend. 
However, deeper and more accurate photometric data as 
well as a treatment of the orientation of the galaxy will 
better constrain these suggestions. 



5. Summary and Conclusions 

In this paper we report observations in the /, J and K s 
band covering about 40' x 30' of the Draco galaxy. The 
photometry has been calibrated to the DENIS catalogue 
in the I band and to the 2MASS catalogue in the J 
and K s bands. The cross-identification with catalogues 
of confirmed members, non-members, carbon and RGB 
stars has helped us to to evaluate the effect of the fore- 
ground contribution which has also been discussed using 
2MASS stars in a comparable field 0°.4 North of the cen- 
tre. Many RGB stars have been detected as well as can- 
didate AGB stars. The TRGB location in the /-band has 
been used to derive a distance modulus for the galaxy of 
(m - M) = 19.49 ± 0.06(stai) ± O.lb(sys) in agreement 
with values in the literature and the most recent deter- 
mination from RR Lyrae stars (Bonanos et al. 12004(1 . The 
surface distribution of all detected sources traces a smooth 
and elliptical structure with e = 0°.27 and PA sw 90° 
also in very good agreement with previous measurements 



(Hodge 1964, Odenkirchen et al. .2001), and suggests a 
variation of the inclination with radius. 

The study of the distribution of stars in the colour- 
magnitude diagram revealed a population with an aver- 
age [Fe/H] of —1.95 dex (or —1.34 dcx) with a spread 
of 1.26 dex (or at least 0.72 dex) using the Valenti et 
al. (2004) photometric indexes (or by comparison with 
the fiducial RGB tracks of a few Galactic clusters). The 
Girardi et al. I|2000) isochrones of two populations with 
[Fe/H]= —1.3 and [Fe/H] — 1.7 better represent the colour- 
magnitude diagrams of sources located in three concentric 
regions. They show a non-negligible intermediate-age pop- 
ulation even though most of the stars are as old as the age 
of the Universe in agreement with Aparicio et al. ((2001) 
who derived that some 75 — 90% of the stellar population 
of Draco is older than 10 Gyr. 
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